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ABSTRACT 


Along with the development of powerful microprocessors and 
microcontrollers, the applications of the model predictive controller, which 
requires high computational cost, to fast dynamical systems such as power 
converters and electric drives have become a tendency recently. In this paper, 
two solutions are offered to quickly develop the finite set predictive current 


control for induction motor fed by 3-level H-Bridge cascaded inverter. First, 

the field programmable gate array (FPGA) with capability of parallel 
Keywords: computation is employed to minimize the computational time. Second, the 
hardware in the loop (HIL) co-simulation is used to quickly verify the 
developed control algorithm without burden of time on hardware design 
since the motor and the power switches are emulated on a real-time platform 
with high-fidelity mathematical models. The implementation procedure and 
HIL co-simulation results of the developed control algorithm shows the 
effectiveness of the proposed solution. 
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1. INTRODUCTION 

Nowadays, most commercial inverters for AC machines are installed with field-oriented-control 
(FOC). To further improve the dynamic response of the FOC, new control strategies have been studied for 
the current control loop such as synchronous vector control, state feedback control, deadbeat control, neural 
network and fuzzy control [1]. Among them, model predictive control (MPC) has been considered as a 
powerful alternative control method for power converters and electrical drives [2-8]. The MPC, which was 
early developed in the 1960s, is an application of the optimal control theory [2], in which the system model is 
used to predict the future behavior of the system states in a predefined time horizon and then obtain the 
optimal control action which minimizes a given cost function. Despite of advantages such as intuitive 
concept, quick dynamic response, ability to handle constrained linear and nonliner multivariable dynamic 
systems [3], the application of MPC is restricted in the field of process control due to its relatively high 
computational cost. In recent years, with the development of of high-speed microcontroller and FPGA, the 
high computational cost problem of the MPC can be solved following that the application of the MPC has 
been expanded to fast dynamical systems such as power electronic converters and electrical drives. Due to 
the fact that the power electronic converters only have finite switching states, finite control set model 
predictive control (FCS-MPC) has been developed to select the switching state for the converters directly 
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instead of using a modulator which is also complicated, especially in multi level inverters. The FCS-MPC has 
been successfully applied to the current control in three-phase inverters [9-14], matrix converter [15]. And 
induction machine [16-18]. An improved FCS-MPC is also proposed to solve the variable switching problem 
ininduction motor (IM) fed by 2-level voltage source inverter (VSI) [19]. 

Multilevel converters is widely used for medium-voltage (HV) and high power applications due to 
its remarkable advantages over the conventional two-level VSI such as: lower voltage stress, lower rate of 
voltage change (dv/dt), reduced total harmonic distortion (THD) and lower switching frequency which results 
in improved switching loss [9,20]. Among the multilevel inverter topologies, the CHB is the most successful 
configuration and already commercialized because of its modularity [21-23]. Typically, a motor control 
system consists of a conventional controller such as proportional-integral-differential (PID), state feedback, 
etc., and a pulse width modulation (PWM). For multilevel inverters, the higher the level, the more complex 
the PWM algorithm is, especially when other inherent problems of the multilevel topologies such as 
capacitor voltage balance, fault tolerant ability are considered. Hence, this research focuses on the application 
of FCS-MPC to IM fed by multilevel cascaded H-bridge (CHB) inverters. Instead of using the conventional 
PWM, the voltage vector which minimizes the predictive current tracking error is selected directly in every 
consecutive sampling cycle. Without PWM algorithm and by using FPGA [24-26] as the controller, the FCS- 
MPC algorithm can easily be sovle in very short time following that the response of the inner current is 
significantly improved in comparison with the conventional PID controller. Besides, multiple objectives can 
also be achieved by appropriately choosing the cost function. To avoid wasting time on hardware design and 
the risky experiments with high power systems, HIL platform manufactured by Typhoon is employed to 
quickly verify the developed control algorithm. By using the Typhoon HIL, the behavior of the motor and the 
CHB inverter are precisely evaluated in real-time by high-fidelity mathematical models. Realtime 
simulations with various scenarios are conducted and the results show that the employed FCS-MPC applied 
to CHB-fed IM drive can achieve good performance not only in terms of tracking accuracy but also 
dynamic response. 


2. SYSTEM DESCRIPTION 
2.1. Overview of CHB inverter 

The typical configuration of a three-phase multilevel CHB is shown in Figure 1 (a). The 
fundamental component of the CHB is a single-phase H-brigde inverter which is normally called a power 
cell. Each phase of the CHB consists of several cells connected in series. 





Figure 1. (a) Structure of a CHB inverter; (b) Topology of a cell 


Assume that each phase consists of n cells connected in series and each cell is fed with an 
independent dc voltage source Vac. The power switches S,-S4 in cell i of a phase can be controlled to generate 
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three voltage levels {+Vac, 0, -Wac} corresponding to cell state S,; ={1,0, -1} where x represents phase a, b or 
c. The voltage level of a phase is defined as the sum of all cell states in that phase as following: 


5, = Sy (1) 


The number of possible voltage levels for one phase 1s: 
m=2n+1 (2) 


The output voltage of each phase is calculated by 


Van = Vo Ss Vegy — Vi Siva Vi S (3) 
i=0 i=0 


i=0 


2.2. Space voltage vector of 3-level CHB inverter 
A voltage vector is formed by a combination of switching states. The total combinations of the CHB 
inverter is as follow: 


K =(2n+1) (4) 
K, =12n° +6n+1 (5) 


In multilevel CHB converters, same voltage vector can be implemented by several combination of 
switching states. Hence, the number of voltage vectors are normally less than the total switching states as 
following [9]. In this research, a three phase three level CHB inverter is used with 27 possible voltage vectors 
in static aß coordinate. By eliminating the switching states which generate high common-mode voltage, the 
number of employed voltage vectors are reduced to 19 as shown in Figure 2. These voltage vetor are 
employed in the FCS-MPC in the next section. 





Figure 2. Space voltage vector of a 3-level CHB inverter 


2.3. Modeling of the induction motor 
In stationary frame aß, the behavior of the IM can be represented by the following equations: 


S S 


d¥ d¥ 
u = Ri +=; 0= Ri + jz, oY (6) 
dt dt 
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Y’ = Li+ L,i; Y? = L i + Li’ (7) 


mr? 


In which R,,L,R,,L,are resistance and inductance of the stator and rotor, respectively. 
The current and flux of the stator and rotor are expressed in vector form as follow: 


u =u Ug) st, =[i,, il sé =, ip) =l¥.. Viel E E Vel (8) 


Defining new variables as: 








1 L L L 
Paar E S a a e (9) 
L. C, R, An 


A fundamental manipulation on (1) and (2) results in: 


s 


di , 
== At +Bu +PP’, 











dt (10) 
dF: 
TTA, Yi, + Bai, 
t 
with 
lee 0 2a 0 1-ô (l-o)@ wea z i 0 
ÔT, ôT, ôL, ôT, ô T, T, 
i iia Lip ie a eet a 
0 5) i 2 ee ee g we ioe 
ÔT, OT, SL, ô ôT, T, T, 


For MPC control design, the continuous-time model (10) needs to be transformed into discrete-time 
model with sampling period T, by using forward-Euler method as following: 


s — 5 s s 
lokai =a Di, , tru, TDP a 


. (11) 
Tn = D, Pink u A 
in which 
® =I, +T,A;®, =1,+7,A,;1, =T,B; 
1 0 
r, =T,B,; D, =T,P; 1, = 
O 1 
Based on (11), the predictive current in N step ahead can be computed by: 
a lr, 0 . Ofun] TO D, 0 a 0 T 
be nao Or, r; aa 0 u; k+2 + o os @D, D, rire 0 Yin k+1 
= ‘ i , : . lok : : A : 
tins `T, `T sar r; Wein p` 0” 'D, ©” D,a a Dani TA 
(12) 


3. CONTROL DESIGN 
3.1. FCS-MPC design 

The block diagram of the IM control system is shown in Figure 3(a) which conists of two outer 
loops including a speed regulator, a flux regulator and an inner current loop using FCS-MPC. The two outer 
loops designed in synchronous reference frame dg generate the reference stator currents. Then, these 
reference currents are transformed into af coordinate and used for the FCS-MPC design. Based on the IM 
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model, the optimal voltage vector is chosen to minimize a cost function in which one or several objectives 
may be considered such as: minimal tracking error, capacitor volage balance, minimal switching loss, etc. 
The detail FCS-MPC algorithm which is activated in every consecutive sampling period is shown in 
Figure 3(b), in details. 


Known _ from outer loop 


Calculate reference current i}; (k) 


Measure i,,.(k), olk) 






Calculate stator current i,, (x) 
Flux rotor ®,, (x), angle @.(x) 





4 
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+y controller dq 





























wa) | l POP 
Pralk)*| + Flux [is] / ap | L wi Apply optimal 
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| 
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Waiting for next sampling instant 
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Figure 3. FCS control strategy, (a) IM control scheme, (b) Predictive current control algorithm 


3.2. Cost function selection 
In this research, the cost function is chosen based on the predictive current error in aß coordinate of 
the induction motor. 


3.2.1. Formulation of the current error 

The error between the predictive current and its reference can be expressed as an absoltute value, 
square value or integral value [3]. If only the present tracking error is used by the cost function, absolute 
error and square error give same performance. When the present and past values of the tracking error are 
considered, the square of error gives better tracking performance than the absolute error while the integral of 
error gives the best performance. However, a cost function with integral of tracking error also requires higher 
computational cost. 
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3.2.2. Delay compensation 

In the ideal case, the time needed for calculation is negligible. Assume that the currents are 
measured at time instance tk, the optimal voltage vector that minimizes the error at time instance tk+ı 1S 
chosen and applied immediately at tk. Therefore, the load current tracks the predicted current at tk+. 
However, if the calculation time is significant compared to the sampling time, there is a delay between the 
measured current insant and the applied new voltage vector instant. During the delay time, the previous 
voltage vector is still applied, which makes the load current move away from the reference and increases the 
current ripple. A simple solution to compensate this delay is to take the calculation time into account and 
apply the selected voltage vector after the next sampling instant [3, 14, 16]. The load current reaches the 
predicted value at tķ+2. 


3.2.3. Prediction of future references 
In general, the future references are unknown and needed to be estimated by using a second-order 
extrapolation [10]. However, for sufficiently small sampling time, the future reference value at time tk+2 1S 


assumed to be approximately equal to the present reference value at time t, [3]: Li (k) = ig (k+2). Based on the 


aforementioned analysis, the cost function is chosen as: 


J= 





i, (k)-i, (k +2)|+ 





is (k) -ip (k +2) (13) 


ig (k) : Reference stator current of the motor in aß coordinate at time instance tk. 


i,,(k +2): Predictive stator current of the motor in aß coordinate at time instance tk+2. 


4. FPGA IMPLEMENTATION 
4.1. Functional internal circuit (IC) design. 

A complicated algorithm can be divided into many smaller calculation steps and an internal circuit 
needs to be designed and programed functionaly to perform each step. An IC includes two fundamental 
blocks: a finite state machine (FSM) and a processing unit (PU) as shown in Figure 4(a). The FSM block 
contains all finite states of the IC and control signals to processing unit. Meanly or Moore method are usually 
used to design the FSM with several typical signals: clk, rst, init and done. In which, Clk and rst are operating 
clock and reset signal for the circuit, respectively. Init and done signal is set to ‘1’ for only one clock period 
when the IC starts and finishes its operation. The done signal is connected to the init signal of the other 
circuit so that sequential calculation can be implemented. Since an internal circuit only operates in a fixed 
small time and is inactive most of the time, this design avoids the propagation of unespected glitches and 
reduces the FPGA power consumption. The PU block contains operators +, -, x, / and takes Data_in as input 
data to calculate and then send the results to the Data_out under the control of the FSM. The structure of the 
PU is designed by the pipelined structure as shown in Figure 4(b) to synchronize the data and signals. c. 
Predictive current control implementation based on FPGA platform. 

The flow chart shown in Figure 5(a) illustrates how to implement the FCS-MPC for current loop of 
the IM on a FPGA platform. To solve the algorithm described in the previous section, there are nine steps 
corresponding to nine internal circuits. The name of Step! to Step9 are: ADC_read, abc_to_ af, Dq to oP, 
Is_to_flux, Pre_model, J_calc circuit, Find_mindJ, af to dq and Flux_ model, respectively. Since Step2 
requires the data calculated in Step/ while Step3 requires the data from Step2, the operation of circuits 1,2,3 
is sequential. Similarly, the circuits 3,5,6,7 and 8,9 must run sequentially. Meanwhile, Step5 requires the data 
from both Step3 and Step4, so that circuit 3 and 4 must operate in parallel. Because circuit 4 needs more 
calculation time than circuit 3, circuit 5 only take the input data when circuit 4 ends its operation, which 
means the done signal of circuit 4 is the init signal of circuit 5. On the other hand, Step9 only requires data 
from Step which means circuit 8 and 9 operate in parallel with circuit 3,5,6,7. Finally, the execution time of 
the FCS-MPC is shown in Figure 5b. As can be observed, several steps can be carried out in parallel which 
reduces the computational time of the implemented algorithm. The parallel computation ability of the FPGA, 
which is not available in conventional DSP/microcontroller, allows the MPC to be applied to many other 
fields without having to worry about the computational time. 
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Figure 4. Functional, (a) IC structure, (b) Pipelined structure 
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Figure 5. Flowchart, (a) FPGA implementation (b) IC’s execution time 
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4.2. HIL-FPGA platform 

To quickly verify the FCS-MPC for IM without risky experiments, real-time simulation based on 
HIL platform are conducted. The virtual IM motor with parameters provided in Table 1 and the 3-level CHB 
inverter are simulted by Typhoon HIL 402 device. The control algorithm is implemented by a FPGA card 
named ZYBO-27. The block diagram and the experimental devices of the real-time simulation system are 
shown in Figure 6 and Figure 7, respectively. 


3-level 
converter 


CHB 













Predictive (Sade 
current at 
instant k+2 { 








Lisl 
Pi ap (x) iap(k) 


Flux model 


Figure 6. Block diagram of the HIL-based simulation system 


Table 1. Parameters of the IM 


Symbol Quantity Value 
Prate Rate power 2.2kW 
Ts Rate torque 7.3Nm 
Liate Rate phase current 4.7A 

V iate Rate phase voltage 400V 
frate Rate frequency 50Hz 
Niate Rate speed 2880 rpm 
R, Stator resistance 1.99 Q 
R, Rotor resistance 1.84 Q 
La Mutual inductance 0.37 mH 








Figure 7. Experiment devices 


5. RESULTS AND DISCUSSION 

In this section, various real-time simulations are carried out to verify the developed FCS-MPC 
applied to IM control. First, the transient response of the control system with the load disturbance is tested. 
At time instance t=/{0.05s, 0.5s, 0.75s}, the load torque is suddenly changed with corresponding values 
T,={0, 0.5T,ate , Trae } while the reference speed is fixed at w,.=300(rad s71). As can be observed in 
Figure 8, the electromagnetic torque generated by the IM quickly tracks the load torque to maintain the rotor 
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speed at its desired value. From time instance t=/.5s, speed reveral of the IM with constant load is 
investigated. The reference speed is changed at time instance t=/.5s and t=4.3s with corresponding values 
Nref=-Nrate and N,er-O0. In this case, the rotor speed also quickly trackes its reference with negligible tracking 
error, 1.e., about 0.2 rad/s. In the second scenario, the ability to generate electromagnetic torque at standstill 
of the IM control system is clarified. With employed FCS-MPC for the current loop, the generated 
electromagnetic torque quickly tracks the load torque at zero speed in just 3ms as shown in Figure 9. This 
quick response is sufficient for most practical applications and comparable to the well-known direct 
torque control (DTC). 
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Figure 8. Electromagnetic torque, (a) Time-varying speed motor and load torque, (b) Zoom-in result 
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Figure 9. Generated torque, (a) Load torque and load current at standstill, (b) Torque step response 


6. CONCLUSION 

In this paper, an advanced solution to quickly develop the FCS-MPC for IM fed by multilevel CHB 
inverter is discussed. First, the behavior of the stator current in future is predicted with each of possible 
voltage vector based on the explicit discrete-time mathematical model of the IM. The optimal voltage vector 
which fulfill the objective of a predefined cost function is selected and applied to the CHB inverter. Second, 
the FPGA with parallel calculation ability is used to minimize the computational time of the complicated 
control algorithm like MPC. Finally, to eliminate the burden of time on hardware design as well as to avoid 
potential risks with high power experimental systems, real-time simulations with high-fidelity mathematical 
models implemented by HIL platform are conducted to verify the control algorithm. Achieved real-time 
simulations show that the developed FCS-MPC gives good performance not only in steady-state but also in 
transient-state. 
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